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py/j ! Abstract 

■ It is well known that the distributions of the interiors of the typical cell of a 
1^ . Poisson line tessellation and a STIT tessellation with the same parameters coincide. 

! In this paper, differences in the arrangement of the cells in these two tessellation 

models are investigated. In particular, characteristics of the set of cells neighbouring 
the typical cell are studied. Furthermore, the pair-correlation function and several 
mark correlation functions of the point processes of cell centres are estimated and 
^ [ compared. 
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1 Introduction 

^ ! Random tessellations form a versatile class of models in stochastic geometry and spatial 
I statistics. They are applied to such diverse fields as the modelling of cellular materials, road 
systems or animal territories. Having the variety of applications in mind, it is of particular 
interest to compare different tessellation models and to provide criteria for good model 
choices. In this paper we are interested in a comparison of two basic tessellation models: 
Poisson hyperplane tessellations and STIT tessellations. Poisson hyperplane tessellations 
are obtained by a subdivision of the Euclidean space M'^, d > 2, into polytopes by random 
hyperplanes following a Poisson law [9] . Iteration stable random tessellations (called STIT 
tessellations for short) were introduced more recently in [6] and can be regarded as outcome 
of a random process of repeated cell division which makes them attractive for particular 
applications. 

A comparison of these two tessellation models is of special interest, because the distri- 
butions of the interiors of their typical cells coincide. Hence, straightforward characteristics 
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such as the cell size distribution or contact distributions cannot be used to distinguish these 
models. For this reason, characteristics describing the arrangement of cells in both mod- 
els are of interest. In [10], [15] the second-order moment measure of the random length 
measure concentrated on the edge system of planar STIT tessellation was investigated. 
It turned out that this measure coincides with that of a suitable Boolean model of hne 
segments having a particular length-direction distribution. Moreover, this second-order 
characteristic was shown to be rather close to but different from that of a Poisson line 
tessellation with the same parameters. 

In order to analyze differences between the STIT tessellations and other well established 
tessellation models, second-order quantities of the systems of nodes, edges, and facets of 
two- and three-dimensional Poisson hyperplane, STIT and Poisson- Voronoi tessellations 
were compared in [8] . It was shown that STIT tessellations show an intermediate behaviour 
between the two other models. 

In this paper, we study the arrangement of cells in planar STIT and Poisson line 
tessellations using characteristics of the cells neighbouring the typical cell. Furthermore, 
second-order characteristics of the point process of cell centres are considered. Since most 
of these characteristics are (currently) intractable to analytical investigation, they will be 
studied using Monte-Carlo simulations. 

We start by introducing the required notions and notation. In Section 3, general rela- 
tions for the neighbours of the typical cell in a random tessellation are derived. Numerical 
values for the means of certain characteristics of the neighbours of the typical cell in Pois- 
son line and STIT tessellations are obtained analytically or by simulation. In Section 4, we 
study the point processes of cell centres in both models by means of their pair- and several 
mark- correlation functions. The paper closes with a discussion of our results in Section 5. 

2 Tessellations and tessellation models 

2.1 Tessellations 

A collection M = {Pj)j>i of polytopes Pj, called cells in our context, is a tessellation of 
if 

(i) M is locally finite, i.e. any bounded subset of M!^ has non-empty intersection with 
only finitely many elements of M, 

(ii) the cells of M have pairwise disjoint interiors, i.e. for i ^ j, int(Pi) fl int(Pj) = 0, 

(ni) the cells of M cover the whole space R'^, i.e. Up^eM ~ 

We denote by Jv[ the space of tessellations in M'^, which can be equipped with a suitable 
cr-field 9Jt, cf. [9], [12]. Thus, a random tessellation is just a random variable taking 
values in the measurable space (A^,OJl). 

A random tessellation is called stationary if its distribution is invariant under transla- 
tions of Moreover, it is called isotropic if its distribution is invariant under rotations 
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around the origin. Prom now on we will assume that any of our tessellations is stationary, 
but not necessarily isotropic. We will indicate whenever we need this additional assump- 
tion. 

2.2 Poisson hyperplane tessellations 

A hyperplane H in R*^ can be parametrised by a pair {p,v) G M x S"^"^, with S"f~^ being 
the upper {d — l)-dimensional unit half-sphere. Geometrically, p is the signed distance 
of H to the origin and v denotes the unique normal vector of H belonging to S"^^. Let 
TZ he a probability measure on S^^, satisfying span(supp(7^)) = W^. Denote by A the 
Lebesgue measure on R and consider a Poisson point process on R x with intensity 
measure 7A <S> TZ, where < 7 < cxo is a constant. Identifying the points of this pro- 
cess with hyperplanes, a Poisson hyperplane process in is obtained. The collection of 
random polyhedra generated by a Poisson hyperplane process is usually called a Poisson 
hyperplane tessellation, cf. [9], [12]. Below, our attention will often be restricted to the 
planar case d — 2, where such a random tessellation is called a Poisson line tessellation. 
If 7Z is the uniform distribution on S*^, the Poisson line tessellation is also isotropic. A 
realisation of a stationary and isotropic Poisson line tessellation is shown in Figure 1. 

2.3 STIT Tessellations 

STIT tessellations formally arise as limits of rescaled iterations (also called nestings) of 
stationary random tessellations in the plane and were formally introduced in [6]. Locally, 
this is within bounded windows W C M*^, STIT tessellations may be interpreted as outcome 
of a spatio-temporal process of subsequent cell division, which is roughly described as 
follows. Fix i > 0, a distribution TZ on S"^^ satisfying the same assumption as in the 
previous subsection and a polytopc W C W^. Now, a random life-time is assigned to W 
which is exponentially distributed with parameter related to the integral-geometric mean 
7^-width of W . After this random life time has expired, a random hyperplane H with 
normal direction drawn from TZ is thrown onto W . The {d — 1) -dimensional polytope 
HnW {di, line segment in the planar case d — 2), splits W into two new cells W+ and 
Now the cell splitting process starts anew and develops independently in both cells W+ 
and i.e. W+ and W- are provided with independent and exponentially distributed 
random life times with parameter related to their mean widths. The fact that the facets 
born during this construction are always chopped off by the boundary of their mother-cell 
yields the characteristic non face-to-face situation (see below) . 

It can be shown that for any time t > 0, the distribution of the tessellation M constructed 
inside W is independent of W . Hence, there exists a whole space tessellation extending 
M. This tessellation is stable with respect to iteration, a property that explains the 
abbreviation STIT. A simulation of a stationary planar STIT tessellation is shown in 
Figure 1. 

Below, the characteristics of a STIT tessellation are marked by a superscript S to 
distinguish them from those of a Poisson hyperplane tessellation marked by superscript P. 
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Figure 1: Realisations of an isotropic Poisson line tessellation (left) and an isotropic STIT 
tessellation (right). 



2.4 Properties of Poisson hyperplane and STIT tessellations 

Let us first mention geometric interpretations of the model parameters 7 and 7^ of a Poisson 
hyperplane tessellation, as well as of the parameters t and 7^ of a STIT tessellation. The 
constant 7 is the surface density (this is the mean [d— l)-dimensional surface area per unit 
d-volume) of the hyperplane tessellation, i.e. in the planar case 7 = L^' and j = Sy for 
(i > 3, as it is typically denoted in stereology. In the STIT tessellation, the construction 
time t also equals the surface density, i.e. t = ii d — 2 and t = Sy for d > 3. In 
both models, the distribution TZ is the distribution of the normal direction in a randomly 
chosen facet point of the tessellation (note that this is not the typical facet, but it may be 
interpreted as the area weighted typical facet). 

Poisson hyperplane tessellations and STIT tessellations share another important property. 
In [6] it was shown that the distributions of the interiors of their respective typical cells 
coincide, whenever we chose 7 = t and the same directional distribution TZ. Recall that the 
typical cell of a stationary random tessellation is, roughly speaking, a random cell uniformly 
chosen within a large observation window. The fact that the typical cell distributions 
coincide may be interpreted by saying that Poisson hyperplane and STIT tessellations can 
be built by using the same pieces of a puzzle. However, the spatial arrangement of these 
pieces is quite different in both cases. For instance, any vertex of a Poisson hyperplane 
tessellation has 2d pairwise coUinear outgoing edges, while vertices of STIT tessellations 
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have exactly d+1 adjacent edges and exactly two of them are coUinear. It is our purpose 
to further analyze the difference of the spatial arrangement of cells of the two different 
tessellation models in the particularly interesting planar case from different perspectives. 

2.5 Typical and weighted i-plates 

The faces of a convex polytope P are the intersections of P with its supporting hyperplanes. 
We call a face of dimension i, i E {0, . . . , d}, an i-face of P. In the special case i = we use 
the term corner, for i = 1 the term side, faces of dimension d—1 are called the facets of 
P and P itself is its only d-ia,ce. We write J^i{P) for the set of z-faces of a polytope P and 
jFi{M) — Upg^ J^j(P) for the set of i-faces of all cells of the tessellation M. Furthermore, 
let F{y) be the intersection of all cells of the tessellation containing the point y. Then F[y) 
is a finite intersection of d-polytopes and, since it is non-empty, F[y) is an i-dimensional 
polytope for some z G {0, . . . , d}. Therefore, we may introduce 

Mi^{F{y) : d\uvF{y) ^ i,y eW}, i^O,...,d, 

the set of i-plates of the tessellation M. Then, an i-ia.ce H G J^i{P) of a cell P E M 
is the union of all those i-plates of the tessellation contained in H. A tessellation M is 
called face-to-face if the faces of the cells and the plates of the tessellation coincide, i.e. 
if Mi = Ti{M) for alH = 0, . . . , d. 

The system Mj forms a stationary particle process of i-dimcnsional polytopes in the 
usual sense of stochastic geometry. The intensity of this process is denoted by Aj. For 
the models studied in this paper we will always have < Aj < oo for z = 0, . . . , d. In 
the planar case, the elements of Mq, Mi, and M2 are the vertices, edges, and cells of 
the tessellation, respectively. Note that the number of vertices or edges of a cell is not 
necessarily the same as the number of its corners or sides, respectively. 

Using the Palm distribution of the particle process Mj, we can define the typical i- 
plate of a random tessellation, see [9], [12] for the technical details. For any of the classes 
Mj, i — 0, . . . ,d, introduced above we will denote by Tj the typical object of the respective 
class. Informally speaking, Tj is what we get by equiprobably choosing an object from Mj 
in a large bounded region of the tessellation. The expectation with respect to Tj will be 
denoted by Ej. 

Besides typical i-plates, one can also consider weighted z-plates where an i-plate 
is selected according to an individual weight. For instance, let nj{x) be the number of 

j-dimensional plates of x E Mi {0 < j < i < d) and let Vj{x) denote the sum of their 
j-volumes. Attach to any x G Mi the weight rijlx) or Vj{x). By picking a random x E Mi 
according to these weights, we obtain Tj[nj] and Tj[vj], the nj- or the t^j- weighted typical 
i-plate of M, respectively. For example rrf[no] is the vertex number weighted typical cell, 
whereas Ti[vi] is the length weighted typical edge. Expectation with respect to Ti[nj] and 
Ti[vj] will be denoted by Ej^„^. and Ej^^,^.. 

We let Nj{Ti) = Ejnj(x) be the mean number of j-plates adjacent to Tj and Vj{Ti) = 
KiVj{x) be the mean j- volume of all j-plates of Tj. In the same manner, we can define 
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the weighted mean values Nj{Ti[nk\) — Ei^„^nj(a;) and Vj{Ti[vk\) — Ej^^^Vj(a;). Moreover, 
we let Nij be the mean number of a; G Mj adjacent to Tj. For example A^^.o is the mean 
number of corners of the typical cell or A^o,d is the mean number of cells adjacent to the 
typical tessellation corner. 

In this paper we will adopt the following notational convention: In formulas like Ej/(x), 
object X is the random i-dimensional polytope the expectation Ej refers to. 

3 Neighbours of the typical cell 
3.1 General formulae 

Let M be a stationary random tessellation in R'^. We say that two i-plates of M are 
neighbours, if they share a common (^ — l)-plate. For x G Mi we denote by J\f(x) = {y E 
Mi : y n a; G Mj_i} the collection of its neighbours. In particular, we will be interested in 
characteristics of the cells neighbouring the typical cell of M. 

We start by noting that A^d-i(^d) is the same as the number of neighbours of the 
typical cell, i.e. Na-i{Ta) is the number of elements in J\f(Ta). It is well known (see 
for instance [9], [12]) that Nd-i{T[) = 2'^, in particular A^i(r2^) = Nq{T[) = 4 and 
that Ni{T2) = Nq{T2) = 6, i.e. the typical cell of a Poisson line tessellation has four 
neighbours in the mean, whereas the typical planar STIT cell has six. It appears that 
except N2{T^) = 14, the values iVrf_i(T/) are unknown. Moreover, we let iVf " and V^'"" 
be defined as 

i.e. Nf''^ is the mean sum of the number of i-plates of cells neighbouring the typical cell, 
whereas V^''^ is the mean sum of i- volumes of the i-skeletons of the cells neighbouring the 
typical cell (each time counted with multiplicities). 

It can be shown that the mean values Nf''^ and y/''* are related to certain mean values 
of weighted typical cells. 

Proposition 1 For a stationary random tessellation we have 

Nt" = 7Vd_i(r,)7v,(r,K_i]), y/'" = iv,_i(rd)i^,(r,K_i]) 

for < i < d. In particular 

Nf'^ = No{T2)NiiT2[no]), 1/,''' = NoiT^MT^im]) . 

Proof. The first line corresponds to Corollary 5 in [14] and the second line is a special 
case of the first when taking into account that the number of edges of the typical cell is 
the same as the number of its vertices in the planar case. □ 

The next result appears in the special face-to-face situation as Theorem 10.1.1 in [9] and 
is a main tool to establish distributional equalities for random tessellations. 
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Theorem 1 For a stationary random tessellation M, let h{M, x) and f{M, x) be translation- 
invariant m,easurahle plate characteristics defined for x G Mj, i = 1, . . . ,d, i.e., h{M + 
y^x + y) = h{M, y) and f{M + y,x + y)^ f{M, x) for all y eR'^. Then 



E, 



h{M,x) J2 fi^^y) 



= E, 



f{M,x) HM,y) 



(1) 



Proof. In the general case, the assertion is obtained as a combination of [2, Theorem 
3] and [4, Theorem 5.1]. For the special face-to-face situation we cite [9, Theorem 10. !.!].□ 

Proposition 1 tells us that certain mean values for the neighbourhood of the typical cell of 
a tessellation can be expressed in terms of mean values of the {d — l)-dimcnsional plate- 
number weighted typical cell (which is the same as the vertex number weighted typical cell 
in the planar case). However, much more is true: 

Theorem 2 Let M be a stationary random tessellation. Then 



E,; 



ceAf(x) 



:E,iVi_i(x)Ei,„,_J(M,x), 
where f is a plate characteristic as in Theorem 1. 



Proof. Using (1) for x & Mi and = 1 we get 



yeMix) 



E,[/(M,y)#A^(y)] = E,iV,_i(x)E,,„,_ J(M, x). 



where the last equality follows from the definition of the weighted mean values. This com- 
pletes the proof. □ 

The statement of the last theorem may be rephrased as follows: Geometric character- 
istics summed up over the collection of all neighbours of the typical cell of a stationary 
random tessellation are the same as geometric characteristics of the {d — l)-plate number 
weighted typical cell multiphed by the number of neighbours of the typical cell of the tes- 
sellation. However, this should not mistakenly be understood as the statement that an 
equiprobably selected cell from the family of all neighbours of the typical cell has the same 
distribution as the {d — l)-plate number weighted typical cell. This is due to the fact that 



= E, 



nd-i{x) f-f 



ceN{x) 



(2) 
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is not necessarily the same as 



'd,d 



1 



(3) 



3.2 Numerical values 

In this subsection we restrict our attention to the planar case and consider a Poisson line 
tessellation with edge length density La and directional distribution TZ. Since TZ can be 
extended to an even probability measure on the whole unit circle S^, a natural zonoid 11 is 
associated with such a tessellation (see [9]). Its polar body in the classical sense of convex 
geometry is denoted by 11°. The general formulas from [3] allow us to express the mean 

2 2 2 2 / \ 

values N- ' and ' by the area of 11 and 11°. In the isotropic case, we have V2[U.) — 
and 1^2(11°) = j^, which leads to the explicit values given in Table 1. 



Table 1: Mean values of the sums of characteristics of all neighbours of the typical cell of 
a Poisson line tessellation. 

No corresponding exact formulas are available for STIT tessellations. Even if their typ- 
ical cells have the same distribution as the corresponding cells in Poisson hne tessellations, 
this is not true for the vertex-number weighted typical cell. This is due to the crucial 
non face-to-face situation and the fact that STIT cells have additional vertices on their 
boimdaries which are no corners of it. 

Considering geometric cell characteristics summed over all neighbours of the typical cell 
is somewhat unsatisfactory with respect to our aim of comparing Poisson line tessellations 
with STIT tessellations, because the mean number of neighbours of the typical cell in both 
models differs. For this reason it is more natural to average these characteristics and to look 
at the 'typical neighbour' of the tessellation. This is closely related to a discussion initiated 
in the context of Voronoi tessellations in [13], where numerical differences between / 
and f'^''^ given by (2) and (3), respectively, were considered for the planar Poisson Voronoi 

tessellation. It is currently unclear from the theoretical point of view how the values / 

^ 2 2 2 2 

and /'^''^ are related to each other. Here, we are interested in the values of V2 and V^' . 

2 2 

Furthermore, the characteristic Cq' , the sum of the number of corners of the neighbours of 

2 2 

the typical cell, and its mean Cq arc investigated. Note that this characteristic coincides 

2 2 

with A^n' in Poisson line but not for STIT tessellations. 



general isotropic numeric 
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Since the mean value characteristics are analytically intractable in both considered 
models, they have to be studied by simulation. For this purpose, the neighbourhood 
systems of 15,000,000 simulated cells of both the STIT and the Poisson line tessellation 
were evaluated. The cells were obtained by repeated simulation of the tessellations in the 
window [—100, 100]^. To avoid edge effects, the cells included in the statistics were selected 
using a minus samphng edge correction. The estimated values of Cg'^, V^''^, and V^'^ as 
well as the means according to (2) and (3) are summarised in Table 2. Our simulations 
underpin the observation already made in [13] and suggest that / ' > f'^''^. The ratio 
between both values in the considered cases is between 0.94 and 0.98. Furthermore, it 
turns out that the 'typical neighbour' in a STIT tessellation is larger w.r.t. area, perimeter 
and number of corners than that in a Poisson line tessellation. 





STIT 


PLT 


PLT theoretic 


^0 


25.865 


16.935 


16.935 


7^2,2 


4.386 


4.357 




/^2,2 


4.311 


4.234 


4.234 




30.848 


15.459 


15.503 


vT 


5.471 


4.078 






5.141 


3.865 


3.876 




49.729 


28.034 


28.070 




8.617 


7.226 






8.288 


7.008 


7.017 



Table 2: Simulated mean values of the sums and means of characteristics of all neighbours 
of the typical cell. 



4 Pair- and mark-correlation functions 

In this section, the arrangement of cells is studied by means of summary statistics associated 
with the point process of cell centres. 

For a stationary point process X with intensity < A < oo, Ripley's function is 
defined as 

K{r):^\¥},[X{B% r > 0, 

where Eq denotes expectation w.r.t. the reduced Palm distribution of the point process X, 
see [12]. If K{r) is differentiable in r, the pair-correlation function g{r) of X is given 

by 
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where Hd stands for the volume of the d-dimensional unit ball. 

For a stationary and isotropic marked point process X (with non-negative real-valued 
marks) we can also define the mark-correlation function as 

E°''"[m(0)m(r)] 

where E°''' is the expectation with respect to the second-order Palm distribution w.r.t. to 
the origin and a point at distance r from the origin, m(0) and m{r) are the marks of these 
points, and n is the mean mark. Hence, k^nm is a normalised mean product of marks of 
points located at the origin and at a distance r from the origin conditioned on the existence 
of these points. 

Pair-correlation functions of the point process of vertices and of the random length 
measure on the edges of the tessellation were studied recently in [8]. In both cases, the pair- 
correlation functions are monotone decreasing. The functions for the STIT tessellation are 
dominated by the functions for the Poisson line tessellation. In particular, the convergence 
of g{r) towards 1 for r ^ cxo is faster for the STIT tessellation. 

Here, we will study pair- and mark-correlation functions of the point process of cell 
centres. Since we are aiming at a comparison of the arrangement of cells whose shape 
distribution is the same, the centres should be depend only on the cell shape, but not 
on the structure on the cell boundaries. Therefore, we chose the center of gravity as a 
centroid function of the cells. Choosing La = 1 and TZ as the uniform distribution on 
100 realisations of both tessellation models were generated in the square W = [—50,50]^. 
In order to avoid edge effects induced by cells hitting the boundary of W, the estimation 
of the correlation functions was restricted to centres within the square [—30,30]^. We 
estimated the pair-correlation function of the unmarked point process of cell centres as 
well as the mark-correlation functions of the process marked with the area, perimeter and 
number of corners of the cells. For the estimation we used the classical translation edge 
correction estimator implemented in the R-package spatstat [1]. The estimated functions 
are shown in Figure 2. 

The plot for the pair-correlation function indicates a higher degree of clustering of the 
cells in the Poisson line tessellation than in the STIT case. Not surprisingly, the mark- 
correlation functions show that cells with small distance between their center points are 
smaller than average. This phenomenon is more pronounced for the Poisson line than for 
the STIT tessellation. The mark-correlation function for the number of corners differs in 
its structure from the functions for area and perimeter by showing a slight increase above 
the value of 1. The relatively high value for r = stems from the fact that the cells 
have at least three corners. Both tessellation models behave similar with respect to this 
characteristic. 
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Figure 2: Pair correlation function of the point process of cell centers in a STIT (solid) 
and a PLT (dashed) tessellation (top left). Mark correlation functions of the same point 
process marked with the area (top right), perimeter (bottom left), and number of vertices 
(bottom right) of the cells. Note the different scale in the plot for the mark-correlation 
function for the number of vertices. 

5 Discussion 

We have evaluated several parameters, which yield information on the arrangement of 
cells in isotropic planar STIT and Poisson line tessellations. The characteristics of the 
neighbours of the typical cell show that the 'typical neighbour' in a STIT tessellation is 
larger than in a Poisson line tessellation. This is in hne with the observations on the pair- 
and mark-correlation functions for the point processes of cell centers. These indicate a 
higher degree of chistcring in Poisson line than in STIT tessellations. Furthermore, cells 
with nearby centers tend to be smaller in the Poisson line tessellation. These findings 
can be explained by the fact that Poisson line tessellations, unlike STIT tessellations, are 
face-to-face. Hence, the dependence between neighbouring cells is higher in the Poisson 
hne tessellation. 
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We finally would like to point out that similar trends are also expected in higher di- 
mensions. However, it has recently turned out that dimension two is the critical dimension 
for the STIT tessellations (compare with [10] and [11]). This is due to the fact that spa- 
tial dependences in dimension two are rather weak compared with the higher dimensional 
cases and that only in the planar case the global construction of STIT tessellations can 
be split into a negligible warm-up phase and the proper construction phase, which unfolds 
already within a typical STIT environment. Moreover, in dimensions > 3, the number 
of facets and the number of vertices of a cell are no more deterministically related. For 
this reason, mean value formulas and relations become considerably more involved and a 
further weighting procedure of the cells is inevitable. In view of these reasons, we found it 
natural to restrict most of our results to the planar case, where the most explicit results 
are available. 
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